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1. INTRODUCTION 


This report summarizes the research performed at Rice University related to NASA-Ames 
University consortium grant NCC 2-5199 during the two year period 1996 September 1 through 
1998 August 31. The research program, titled Orion Nebula and Planetary Nebulae . involved the 
analysis of Hubble Space Telescope (HST) imagery and spectroscopy of the Orion Nebula and of 
the planetary nebulae NGC 6818 and NGC 6210. In addition, we analyzed infrared spectra of 
the Orion Nebula taken with the Infrared Space Observatory (ISO) The primary collaborators at 
NASA-Ames were Drs. R. H. Rubin, A. G. C. M. Tielens, S. W. J. Colgan, and S. D. Lord (Tielens 
& Lord has since changed institutions). Other collaborators include Drs. R G. Martin (CITA, 
Toronto), G. J. Ferland (U. KY), J. A. Baldwin (CTIO, Chile), J. J. Hester (ASU), D. K. Walter 
(SCSU), and P. Harrington (U. MD). 

In addition to the Principal Investigator, Professor Reginald J. Dufour of the Department of 
Space Physics & Astronomy, the research also involved two students, Mr. Matthew Browning and 
Mr. Brent Buckalew. Mr. Browning will be graduating from Rice in 1999 May with a B.A. degree 
in Physics and Mr. Buckalew continues as a graduate student in our department, having recently 
received a NASA GSRP research fellowship (sponsored by Ames). 

The collaboration was very productive, with two refereed papers already appearing in the 
literature, several others in preparation, numerous meeting presentations and two press releases. 
Some of our research accomplishments are highlighted below. Attached to the report are copies of 
the two major publications. Note that this research continues to date and related extensions of it 
recently has been awarded time with the HST for 1999-2000. 

2. THE ORION NEBULA 

Our research on the Orion Nebula, the nearest and brightest Galactic H II region primarily 
involved UV-optical-IR spectroscopy with HST and ISO. The HST observations involved use of the 
Goddard High Resolution Spectrograph and Faint Object Spectrograph to perform the first study 
of the UV Nil] emission lines at 2143 A &; 2140 A in a Galactic H II region. From measurements 
of these and other UV-optical lines, coupled with detailed photoionization models of the nebula, 
we found that the gas-phase value of N/O is in the range 0.13 to 0.18, a value somewhat higher 
than solar. Moreover, requiring that N+/0 + be the same as the values derived solely from UV and 
optical lines independently, allowed a new value for the rms temperature fluctuation parameter ( t 2 ) 
to be derived in the Orion Nebula. The results were published in the 1998 March 10 issue of The 
Astrophysical Journal (copy attached). They were also discussed in invited presentations by Rubin 
at the Sixth Texas-Mexico Conference on Astrophysics held in Houston during 1997 March and 
by Dufour at the Chemical Evolution from Zero to High Redshift meeting in Garching, Germany, 
during 1999 October. The results have(will) appeared in meeting conference proceedings (e.g., 
Rubin et al. 1998, RMAA , Serie de Conferences, volumen 7, pp. 34-71). 

In addition, we successfully accomplished the first measurements of the infrared He I Bra line 
in the Orion Nebula (or in any H II region for that matter), with the SWS on the ISO spacecraft 
and used it to derive the He^/H* in this nebula. By comparing it to measurements of the nearby 
H I Bra line, we derived a value of He + /H + = 0.085±0.003. This value is similar to values found 
in recent optical studies of Orion, but our technique is less affected by extinction by dust in and 
around the nebula. These results were published in the 1998 July 10 issue of Astrophysical Journal 
Letters (copy attached). 
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3. PLANETARY NEBULAE 

We succesfully obtained imagery observations of two planetary nebulae, NGC 6210 & NGC 
6818, with the WFPC2 CCD camera on HST during 1997. Dufour and undergraduate student 
M. Browning of Rice processed the imagery of NGC 6818 during the summer of 1997 to produce 
calibrated maps of the nebula in the emission of [O III] 4363 A & 5007 A, along with H/3 and Ha, 
at a spatial resolution of about 0.1 arc seconds (about ten times higher than possible with ground- 
based telescopes). The H/? and Ha images were ratioed to produce a map of the extinction by 
dust across NGC 6818, which was then used to correct the map of the observed ratio of the [0 III] 
5007/4363 A line ratio, then derive point-to-point values of the electron temperature, T e , across the 
nebula. This enabled us to study the spatial variations in T e and the temperature fluctuations in 
the nebula. In collaboration with R. Rubin at Ames and P. Harrington at the Univ. of Maryland, 
we used these data to derive a rms temperature fluctuation value for the nebula. We also found that 
the T e variations were irregular across the nebula and did not correlate with surface brightness or 
morphology. The results were presented in a poster (“HST Imagery and Temperature Fluctuations 
for the PN NGC 6818,” by Rubin, Dufour, Browning, and Harrington) at the 1998 January meeting 
of the American Astronomical Society in San Diego (abstract published in BAAS, 30, 896, 1998). 

In addition to the above quantitative mapping of the spatial variations of extinction and 
temperature across NGC 6818, we also constructed color-coded “false-color” images of NGC 6818 
which were subsequently posted in the “Gallery of Planetary Nebulae” at the website of the Space 
Telescope Science Institute (http://www.stsci.edu). During 1998 April Dufour visited NASA-Ames 
to give a colloquium and work with R. Rubin on our processed data on NGC 6818 and newer imagery 
of NGC 6210. Subsequently Rubin and his summer student, C. Ortiz, refined the analysis of NGC 
6210 to produce some rather striking color-coded images of the more complex planetary nebula NGC 
6210, which led to it being called “The Turtle,” and the pictures released by STScI in a press release 
(STScI-PR98-36, http://oposite.stsci.edu/ pubinfo/pr/1998/11/) that was subsequently reported 
by several print and electronic media. Currently, we are working on modeling these two nebula and 
preparing a detailed analysis of their physical conditions and morphology for publication in 1999. 

4. SUMMARY 

This has been a very productive collaboration, as evident by the numerous publications and 
presentations, more of which will be forthcoming beyond the two years of the consortium grant. 
Recently, we have extended our collaborations to include Wolf-Rayet galaxies and the planetary 
nebula, NGC 7009, and expect that new agency funding will be forthcoming in 1999 to assist our 
research efforts in these nebular fields. Noinventions or patents were reported on this grant. 

5. ATTACHMENTS 

(see next pages) 
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ABSTRACT 

Using the Short Wavelength Spectrometer on the Infrared Space Obser\'ator\\ we measured in the Orion Nebula 
several components of the He i Bra transition, the first detections to the best of our knowledge in an H n region, 
and the H i Bra line. A value of He + /H" = 0.085 ± 0.003 is inferred from these data and is typical of previous 
values found for Orion. The IR method to obtain He7H^ is very insensitive to extinction. The measurement of 
the weak He i 2.855 pm (5p } P" -* 4 s \S) flux permits an assessment of T in the He" region. Its flux relative to 
the flux of He i 4.0490 pm (5g 'C-4/'F\ 5 g '0-4/'^" blend) is roughly oc T r (for T f range applicable to 
most H ii regions and planetary nebulae) and depends little on N t .. From our measured ratio, using N r = 10 4 
cm -3 , is slightly less than 10 4 K. We stress that many IR measurements of a hydrogen line are likely to also 
include one or more significant helium line contributions. Attributing all of the flux to the stronger H I line will 
result in underestimates for abundance ratios relative to H~. This has usually not been taken into account in the 
past when deriving abundances of other ions/elements relative to H. Errors in the velocity and line width inferred 
may also occur by failure to consider the blended He component. Because the strongest He I lines will generally 
be on the blue side of their counterpart H I line, interpreting a blended feature as only H i will cause a spurious 
bias to lower velocities. 

Subject headings: H n regions — ISM: abundances — ISM: atoms — ISM: individual (Orion Nebula) 


1. INTRODUCTION 

Gaseous nebulae are laboratories for understanding physical 
processes in all emission-line sources and are probes for stellar, 
galactic, and primordial nucleosynthesis. Most of the helium 
in the universe was produced in the big bang. The reliable 
determination of the primordial He abundance (Y p ) is of fun- 
damental interest and is certainly not a settled issue (e.g., Olive 
et al. 1997). Much of the emphasis on deriving Y p is based on 
measurements of very low metallicity galaxies in which there 
is less buildup of He from galactic nucleosynthesis. While 
knowledge of the helium abundance in a Galactic H ii region 
plays a minimal direct role in evaluating Y p because of the 
relatively high metallicity and the need to extrapolate to zero 
metallicity, there are other reasons why it is important to mea- 
sure an accurate He abundance. For example, here we will be 
testing further the theoretical predictions of the level popula- 
tions for He i (Smits 1996, hereafter S96). Smits’s tables of 
relative line emissivities from a recombination-collision-cas- 
cade model are generally the ones used in deriving He abun- 
dances from He i lines. Because of its brightness and the ability 
to model it with some success, the Orion Nebula is often a 
testbed for applying and verifying various physical processes. 
If some discrepancy from observations is found in Orion, it 
may affect the entire scale of deriving He abundances, including 
those used to obtain Y p . 

If H and He lines are formed by radiative recombination, 
then their relative intensities are ocHe + /H + . However, collisions 
from the metastable 2 3 5 level selectively enhance the intensities 
of some He i lines. As a result, the He abundance inferred 
from such lines will be systematically overestimated unless this 
is properly taken into account. Kingdon & Ferland (1995, 1996) 

1 Based on observations with ISO , an ESA project with instruments funded 
by ESA Member States (especially the PI countries: France, Germany, the 
Netherlands and the UK), and with the participation of ISAS and NASA. 

2 NASA/Ames Research Center. MS 245-6, Moffett Field. CA 94035-1000. 

3 Orion Enterprises. 

4 Rice University, Space Physics and Astronomy, Houston, TX 77251. 

5 IPAC/California Institute of Technology, Pasadena, CA 91 125. 


examined some of these transitions in terms of both collisional 
effects and self-absorption, suggesting a line-by-line modifi- 
cation to the dominant physics included in S96. Also, Oster- 
brock, Tran, & Veilleux (1992) indicated that these effects are 
not large for Orion. 

The main purpose of this Letter is to present our new Infrared 
Space Observatory (ISO) measurements of key H and He lines. 
We use these to infer a He + /H~ abundance ratio. By comparing 
with earlier He + /H* results in Orion, we will consider whether 
we understand the mechanisms for populating the relevant 
atomic levels that govern the IR He i emission addressed. We 
present the basic atomic physics in § 2 and the ISO observations 
and data reduction techniques in § 3. Our results and conclu- 
sions are discussed in § 4. 

2. ATOMIC PHYSICS BACKGROUND 

At high principal quantum number rc, lines produced by an- 
gular momentum splitting, measured by quantum number /, 
blend. At radio wavelengths, recombination lines of H and He 
are adjacent each other in close proximity (e.g., tables in Towle, 
Feldman, & Watson 1996; Lil ley & Palmer 1968). The differ- 
ence in energy between H and He is then essentially determined 
by the different Rydberg constants R u and R He . As n becomes 
smaller, the /-splitting is relatively more important. He i has 
both a singlet and triplet ladder. With Q < l < n — I and the 
selection rule A/ = ±1, there are 4/i — 2 components, where 
n is the lower level quantum number. Thus, for the He i Bra 
transition (n — 5—4), the splitting of the n = 4 and 5 levels 
gives rise to a set of 14 lines. These energy levels and lines 
are shown in Table 1. Several of the individual He components 
can be distinguished from H i Bra with the spectral resolution 
of ISO. 

3. ISO observations 

Spectra were obtained with the Short Wavelength Spectrom- 
eter (SWS) in SWS02 grating mode in 35 bandpasses (some 
overlapped in wavelength coverage) on 1997 October 12 (UT). 


I ">OQ 



L210 


RUBIN ET AL. 


Vol. 501 


TABLE I 


He and H Bra Transitions 




Relative 




Transition 

Intensity 


Flux' 1 ' 

(/im) 

(Upper-Lower) 

Theory* 

Observed* 

<10-* WcnT-) 

2.855020 .... 

... 5p'P'-4s'S 
.. Sp'P’^s'S 
.. Sd 'D-4p 'P‘ 
.. Sp 'P'-id 'D 
.. Sf'F^d 'D 

0.001714 

6 

0.0719 ± 0.0153 

3.330851 .... 

9.264E-4 

No 

3.703571 .... 

0.005595 

No 


4.0063997 ... 

4.287E-4 

No 


4.0377351 ... 

0.02739 

12 

l.l 1 ± 0.134 

4.0409343 ... 

.. Sf'F'-dd 'D 

0.009129 

12 

0.347 ± 0.087 

4.0490143 ... 

.. 5. ? 'G-4/'F’ 
.. 5g 'G-4fF‘ 

0.07561 

12 

4.32 ± 0.1 17 J 

4.0490343 ... 

0.02520 

12 

d 

4.052262 .... 

.. Hi 5*-4’ 


12 

78.0 ± 1.47 

4.0545045 ... 

.. Sd 'D-Af 1 F' 
.. 5</’0-4/'P 

7.48IE-5 

12 

<0.13 

4.0563459 ... 

2.0I2E-4 

12 

<0.13 

4.122730 .... 

.. 5d'D-4p'P' 
.. 5 p’P’-M'D 

0.002223 

13 

<1.01 

4.244067 .... 

0.003122 

No 


4.606601 .... 

.. Ss'S-Ap'P’ 
■ . Ss ’5-4/7 'P‘ 

7.479E-4 

No 


4.694980 

0.001581 

No 



3 With respect to He i 4473 A line for Smits’s TV = 10 4 crrT\ T r = !0 4 K, 
case B. 


b Number entries correspond to the bandpass observed. 

' Uncertainties are I a and are statistical only. When the line is not detected, 
we use 3 a as the upper limit. 

” This flux entry is the sum for these two blended lines (see text). 

SWS02 scans provide the full resolution of the instrument (de 
Graauw et al. 1996). Observations used both the short wave- 
length (SW) and long wavelength (LW) sections of the si- 
multaneous. dual bandpasses available. The total time observed 
was 12,633 s. Most of these data were taken for other purposes; 
here we present data as necessary for this paper. We covered 
H t Bra and eight of the 14 components of He I Bra with three 
separate bandpasses/bands: 6/ IB, 12/ IE, and 13/2 A. 6 For all 
of the Bra lines discussed in this paper, the aperture was a 
rectangle 14* x 20*. The orientation in the sky of the long axis 

* We use the same numbering system used by the ISO Spectral Analysis 
Package (ISAP) to refer to our 35 bandpasses. Data here are from bandpasses 
6. 12. and 13. ISAP is a joint development by the LWS and SWS Instrument 
Teams and Data Centers. Contributing institutes are CESR. IAS IPAC MPE 
RAL. and SRON. 


is at position angle I90°28. The location of the center of the 
aperture is at a, 6 = (05 h 35 m 14!7l, -05°23'4I'.'5) (equinox 
J2000), 18*5 south and 2672 west of 0‘ Ori C. This position 
was chosen to correspond to the center of other apertures used 
to observe with Hubble Space Telescope-FOS-\SW, a circular 
aperture of 0786 diameter and GHRS-ISW, a 1"74 square ap- 
erture (see Fig. I in Rubin et al. 1997). 

3.1. ISO Data Processing 

The data products delivered to us underwent the standard 
pipeline (e.g.. Leech et al. 1997, § 8) using the latest calibration 
files as of 1998 March 17. The detector signals marked in the 
pipeline as “bad data" were then automatically removed within 
ISAP; further removal of dubious data was also performed 
interactively with ISAP. The spectra of the 12 individual de- 
tectors in each observed bandpass were 3 a-clipped averaged 
first over scan direction and over all scans in the specific band- 
pass. We checked whether an application of a detector “flat- 
fielding” using the ISAP “Shift” routine (adjust all detectors 
to the mean signal level) would improve the consistency be- 
tween the 12 detectors. However, this was found not to be the 
case here. Finally, an average over the detectors was performed. 
In all the averaging steps, a rebinning in wavelength was done 
using a bin size ~0.2\/R, where R is the instrumental resolution 
at the particular wavelength (see Sturm 1997 for discussion of 
techniques used here). For bandpasses 6, 12, and 13, we used 
bin sizes of 0.0003, 0.0005, and 0.0005 pm, respectively. 

The ISAP routine Line Fit was used to fit a one- or two- 
component Gaussian to the detectable lines or to set an upper 
limit to the line flux when the line was not detected. Figure 1 
displays the spectrum (bandpass 12) that contains most of the 
components of interest for this paper. Under all physical con- 
ditions expected for Orion, the dominant He i Bra feature is 
the blended triplet and singlet emission lines 5 g } G~*4f i F" 
and 5 g 'G->4f'F" at = 4.0490143 and 4.0490343 pm. 
We use Martin (1987) for energy levels and S96 for line 
strength predictions. According to Smits’s case B recombina- 
tion theory, the line intensity of the triplet is 3 times the in- 
tensity of the singlet (for expected Orion physical conditions). 



4.025 4.03 4,035 4.04 4.045 4.05 4.055 

Wavelength (/im) 


F ‘ C U Portion of SWS bandpass 12 spectrum of the Orion Nebula observed with ISO. All line fits (dashed lines) are made with the final averaeed data (see 
text) but are shown here superposed on all the data points used. The solid line is the sum of the Gaussian fits. Values of X,, for the H t and He t lines of interest 
are marked. The inset provides an expanded intensity scale; points here are the data averaged over scan direction and detectors. 
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If we thus weight the two blended components, the effective 
= 4.0490193 /x m. The difference AX from the H I Bra (at 
4.052262) is 0.0032427 /xm. In fitting the partially blended 
H I and He I lines, we use only one FWHM for both lines 
because the widths are clearly set by the instrumental FWHM. 
The best fit to the He i feature is obtained by fitting it alone 
using the portion that is unaffected by the much stronger H i 
line. The fit to these partially blended lines is shown in Figure 
1, which also provides an expanded scale (inset). The measured 
flux for the blended pair of He i lines in Table 1 is 4.32 x 
10' 19 W cm -2 . The line flux for the H I line is 7.80 x 10" 18 
W cm' 2 , the (instrumental) FWHM is 0.0023 1 /x m, and the AX 
of 0.0032427 /xm was fixed. 

The next strongest He i line is clearly expected to be the 
5 f y F" -* 4 d 3 D component at 4.0377 /xm, which also is meas- 
urable in bandpass 12. The weaker singlet 5f l F° “► 4d { D 
(4.0409 /xm) is also present. We fit each with a single Gaussian 
requiring the line centers to be fixed (to the observed H i center 
minus the known AX’s) and FWHMs to be fixed (to the H- 
line velocity FWHM). See Figure 1 and Table l. The only other 
component of He i Bra that we detect is 5 p y P° -* 4 s 3 S at 
2.8550 /xm in bandpass 6 (Fig. 2). We fit this with a Gaussian, 
fixing only the FWHM = 0.00125 fxn i, provided by ISAP for 
band IB. The central wavelength is 2.8547 ± 0.00008 /xm. We 
can only set an upper limit on the 5 d 1 D 4 p l P° at 4.1227 
ptm component in bandpass 13. As we shall discuss, the 4.1227 
/xm line may be expected to be somewhat brighter than the 
2.8550 /xm line. However, the latter is observed for a much 
longer time and with more sensitive detectors (band IB) com- 
pared with the 4.12 line (band 2 A), resulting in a factor of ~25 
gain in detectability. The remaining two components of He i 
Bra are not detected. These are 5 d 'D —> 4f'F° (4.0545 /xm) 
and 5 d 3 D -+ 4 f y F° (4.0563 /xm), covered by bandpass 12; we 
provide 3 o upper limits for their line flux. 

We use the data in Table 1 without correcting for extinction. 
In the smaller FOS-1SW aperture, we find that the extinction 
correction factors at 2.855 and 4.052 /xm are 1.09 and 1.05, 
respectively. Thus, the maximum uncertainty in ratios of the 
He i components (ranging from 2.855-4.052 /xm) to H i Bra 
due to differential extinction is less than 5%. These extinction 
corrections are based on the observed HSTIFQS data for Bal- 
mer. He i, and [O n] emission lines at position FOS-1SW. This 
set of corrections is consistent with the shape of the extinction 
curve derived from observations of the Orion Trapezium stars 
(Martin et al. 1998). A reddening map (O’Dell, Walter, & Du- 
four 1992) indicates no reason to expect a different conclusion 
for our larger ISO aperture. 

4. DISCUSSION AND CONCLUSIONS 

The measurement of the H i and He i Bra components per- 
mits a new wavelength regime to be used to determine the He V 
H" abundance ratio in Orion. We make the same assumptions 
that are generally made for deriving the He^/TU ratio from 
adjacent radio recombination lines of H i and He i. In particular, 
we assume that local thermodynamic equilibrium (LTE) holds 
for the population distribution in quantum number n. Then the 
departure coefficients {b n -\ alues), which measure departures 
from LTE, are unity. Even if there are departures, the non- LTE 
factors should be the same for both H and He and hence should 
not affect their ratio (e.g.. Brown, Lockman, & Knapp 1978). 

A very useful tool for interpreting our data in terms of a 
He + /H* abundance is S96 and more extensive tables kindly 
provided by D. P. Smits. These list theoretical strengths for the 



Fig. 2. — Portion of SWS bandpass 6 spectrum, showing the He i 2.8550 

line. We plot the Gaussian line tit with a quadratic continuum fit to the 
data, averaged over all four scans and over detectors. 

He I lines for various values of N t . and T r . We use the entries 
for N r = 10 J cm' 3 and T e - 10 4 * K and case B (optically thick 
in the resonance lines), although we shall comment on how 
these change with other values. These are his tabulated values 
that are closest to the expected average values for A, and T r in 
Orion (Baldwin et al. 1991; Rubin et al. 1991). In Table l, we 
list for the He i lines the intensity relative to the 4473 A line, 
as Smits has. The intensity of the blend of triplet and singlet 
(at 4.0490 /xm), which is always the strongest feature in the 
He I Bra multiples is 0.10081. This is 0.65435 of the intensity 
of the sum of all 14 lines in the multiplet. Let us use 
frac 4049 to represent this fraction. With the measured line flux 
of 4.32 x 10“ 19 W cm" 2 , the projected total for the multiplet 
is 6.602 x 10' 19 W cm" 2 . The flux for the H t line is 
7.80 x 10" 18 W cm" 2 ; thus the He i/H i line ratio is 
0.0846 ± 0.0028 (only the flux uncertainties have been prop- 
agated). This may be interpreted as equivalent to the He"/H" 
ratio along the specific observed rectangular column through 
the nebula. It is not in general equal to the He/H abundance 
ratio unless the He 4 and H* ionization structure are the same 
(e.g., see Rubin 1969). 

The value of frac 4049 changes somewhat when different 
Af., T e pairs are considered. Smits has calculated these for N, — 
100, 10 4 , and 10 6 cm" 3 and for T r = 20,000, 10,000, 5000 K 
and lower (continuing by factors of 2). For expected Orion 
conditions, the variation in frac 4049 is slight. Over the full range 
that we have explored, frac 4(>4q varies by only 9% from a low 
of 0.64592 {N t = 10 6 , T e = 10 4 ) to a high of 0.70472 (A, = 
100, T e = 5000). 

Following the above notation, we examine the variation with 
/V and T f of the fractions of the total He i multiplet in some 
of the other key (or measured) lines: frac 40377 , frac 40409 . and 
frac : 855 . The intensity of the 4.0409 line is essentially always 
a factor of 3 weaker than that of 4.0377 over our range of 
interest (case B, S96). Thus, we need comment on only the 
2.855 and 4.0377 lines. For N r = 10 4 cm" 3 , T r = 10 4 K (case 
B), frac 40377 = 0.17779 and frac :8 „ = 0.01 1 126. Over the full 
range that we have explored, frac 40377 varies by a factor of 
1.105 from 0.16435 (A, = 100, T t = 5000) to 0.18163 (A, = 
10 6 , T e = 10 4 ), and frac^ 855 varies by a factor of 1.605 from 
0.0072331 (A, = 100, 7; = 5000) to 0.011606 (A, = 10 ft , 
T e = 10 4 ). 

The variation in both of these fractions is opposite to that 
of frac 4049 in that lower T t and lower A, result in lower values. 
The small variation of frac 40?77 depends on changes in both 
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T t and N t , However, the more substantial variation of frac 2s „ 
depends much more sensitively on changes in T r rather than 
N r . As such, one could envision using the flux ratio of the 2.855 
line to the 4.049 ^m line as a way to measure T r . This ratio is 
predicted to vary by a factor of 1.75 from 0.010264 (/V, = 
100, T r = 5000) to a high of 0.01797 (N, = !0\ T r = 10 4 ). 

Smits (1996) has also calculated the relative emissivities for 
the He I singlets for case A (optically thin in the resonance 
lines). When we use these values instead of case B, there is 
very little difference in any of the frac x values mentioned here; 
for instance, at /V = 10 4 cm"' and 7. = 10 4 K, changes are 
^1.5%. Thus, our further discussion will proceed with case B. 

From Table 1, the observed flux ratio 
7(4.0377 )/F( 4.0490) = 0.257. This ratio is in the range of the 
theoretical tabular values; it is a factor of only 1.057 less than 
expected for N t . = 10 4 , T, = 10 4 K. This close agreement is 
tempered by the uncertainty in the line measurements. Because 
the S/N is much worse for the 4.0377 line than for the 4.0490 
feature, we do not use it in our abundance analysis. 

The ratio of the flux in the weak He t 2.855 pim line to the 
flux in the 4.0490 blend permits an assessment of T c in the 
He + region. The predicted ratio increases monotonically with 
T r \ it is 0.010718. 0.017002, and 0.027048 at 5000, 10,000, 
and 20,000 K, respectively (using N r = 10 4 cm"'). In this T v 
range, which applies to most H n regions and planetary nebulae, 
the flux ratio is roughly oc7^. Thus a measurement of this ratio 
appears to be a promising gauge for planetary nebulae as well. 
From our measured ratio, 0.01665, T t - 9700 K. Although the 
T r value could be more tightly constrained with a finer T e grid 
than provided now by S96, the uncertainties in the 2.855 flux 
make that somewhat academic. 

We compare with other He7H + radio and optical determi- 
nations. Invariably, these will have been made centered at other 
positions in the Orion Nebula and with different size apertures. 
The radio results include 0.091 ± 0.005 (66a line at 22 GHz 
with beam FWHM = 0!7) (Thum, Mezger, & Pankonin 1980). 
Lockman & Brown (1982) found 0.082 ± 0.008 (76a, 10.5 
GHz, 3!2) and 0.088 ± 0.006 (95/3, 14.7 GHz, 2!2) for the 
central lines of sight. They concluded from a review of the 
Orion data available then that there was no correlation of He"/ 
H + with n. Peimbert et al. (1992) found the same result within 
their uncertainties by observing three line pairs with frequen- 
cies —8.6 GHz using the same telescope/receiver (FWHM = 
3!2): 0.0921 ± 0.0025 (9 1 a), 0.0839 ± 0.0052 (114/3), and 
0.1049 ± 0.0080 (I30y). 

There is much optical data for Orion, but observations are 
of considerably more limited spatial coverage compared with 
the radio beams. We summarize the results of two recent stud- 
ies. From their extensive He i data, Osterbrock et al. (1992) 
found He + /H + = 0.0894, and correcting for neutral He, 
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He/H — 0.101. Esteban et al. (1998) did echelle spectropho- 
tometry at two positions. Their final mean values for HeVH* 
are 0.0856 ± 0.0069 and 0.0893 ± 0.0092; their total He/ 
H = 0.0977. We conclude that the He7H + result of 0.085 from 
this study is typical of the values found by others. There is no 
observationally driven imperative here that would indicate a 
departure from the Smits theory. The combined observational/ 
theoretical picture points to the He" volume being smaller than 
the H" volume. In their "blister" model for Orion. Rubin et al. 
(1991) calculated a flux ratio of He i to H 1 radio recombination 
lines -0.083 for the central lines of sight when convolved with 
Gaussian beams up to FWHM - 1\ This resulted when the 
total He/H was 0. 1. The total He/H ratio for Orion is higher 
than the HeVH" ratio, with the difference being attributed to 
some neutral helium in the H 11 region. 

ISO will provide a rich pool of H 1 and He 1 lines in nebulae, 
many with An > 2. Of all those accessible. Bra is probably the 
easiest to use to measure the adjacent pair. It is the brightest 
and has the largest relative separations between their line cen- 
ters. For instance, at Pfa. the wavelength difference between 
H 1 and the strongest expected He 1 component (6// y H" 

5# y G) is 0.003721 pirn. Thus X/AX = 2000. which exceeds the 
range in SWS resolution. R - 1000-1500, and leads to unre- 
solved lines. This contrasts with X/AX - 1250 for the strongest 
Bra lines, while R (-1300-2100) at this wavelength is suffi- 
cient to resolve these components for the most pan (Leech et 
al. 1997, p. 64). Furthermore, the H 1 Pfa intensity should be 
about one-third that of Bra. Already for Pfa. the brightest 
He 1 component shows only as a slight asymmetry on the blue 
wing of the stronger H I line. Indeed, many IR measurements 
of a H r line are likely to also include one or more significant 
He 1 line contributions. By attributing all the flux to the H 1 
line, one will underestimate abundance ratios relative to H". 
This has usually not been taken into account in the past when 
deriving abundances of other ions/elements relative to H. In 
addition, errors in the velocity and line width mav be introduced 
by failure to account for and to correct for the blended He line. 
Because the strongest He 1 lines will generally be on the blue 
side of their counterpart H 1 line, interpreting a blended feature 
as only H 1 will cause a spurious bias to lower velocities. In 
the present case, had the 4.049 He 1 Bra feature remained 
unresolved from the H 1 line, attributed to the H 1 line 
would be 12 km s“' too small. 
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